Under, Over and Critical Damping 

1. Response to Damping 

As we saw, the unforced damped harmonic oscillator has equation 

mx + bx + kx = 0, (1) 
with m > 0, b > and /c > 0. It has characteristic equation 

ms^ + bs + k = 

with characteristic roots 

-b± Vb^ - 4mk 

2m ' ' 

There are three cases depending on the sign of the expression under the 
square root: 

i) b^ < 4mk (this will be underdamping, b is small relative to m and k). 

ii) b^ > Amk (this will be overdamping, b is large relative to m and k). 

iii) = 4mA: (this will be critical damping, b is just between over and 
underdamping. 

Mathematically, the easiest case is overdamping because the roots are 
real. However, most people think of the oscillatory behavior of a damped 
oscillator. Since this is connected to underdamping we start with that case. 

Case (i) Underdamping (non-real complex roots) 

If b^ < Amk then the term under the square root is negative and the 
characteristic roots are not real. In order for b^ < Amk the damping constant 
b must be relatively small. 

First we use the roots (2) to solve equation (1). 
Let co^ = — 4mk\/2m. Then we have 

characteristic roots: — ± icv^. leading to 

complex exponential solutions: e(-2'/2'"+;a;rf)f^ g{-b/2m-ia,a)t _ 

The basic real solutions are e^^*/2m (;os(a;(if) and e^^'/^m sin(c<j(ff). 

The general real solution is found by taking linear combinations of the two 

basic solutions, that is: 

x{t) = cie-^''^"^ cos{LOdt) + cie-^"^"" sin{cOdt) 
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or 

x{t) = e-^'^^^'icicosicodt) + C2sjn{cOdt)) = Ae-^*^^"" cos{wdt - (p). (3) 

Let's analyze this physically. When b = Q the response is a sinusoid. 
Damping is a frictional force, so it generates heat and dissipates energy. 
When the damping constant h is small we would expect the system to still 
oscillate, but with decreasing amplitude as its energy is converted to heat. 
Over time it should come to rest at equilibriim\. This is exactly what we 
see in (3). The factor cos{(x>^t — ^) shows the oscillation. The exponential 
factor e^''f/2m j-^^g ^ negative exponent and therefore gives the decaying 
amplitude. As f — )■ oo, the exponential goes asymptotically to 0, so x{t) 
also goes as5anpotically to its equilibrium position x = 0. 

We call cvd the damped angular (or circular) frequency of the system. 
This is sometimes called a pseudo-frequency ofx{t). We need to be careful 
to call it a pseudo-frequency because x{t) is not periodic and only periodic 
functions have a frequency. Nonetheless, x{t) does oscillate, crossing x = 
twice each pseudo-period. 

Example 1. Show that the system x + Ix + 3x = is underdamped, find its 
damped angular frequency and graph the solution with initial conditions 

x{0) = 1, j(0) =0. 

Solution. Characteristic equation: + s + 3 = 0. 
Characteristic roots: — 1 /2 ± / \/TT/ 2. 
Basic real solutions: e~^^^ cos{Vlit/2), e~*^^sm{VVi.t/2). 
General solution: 

x{t) = e-*^^{cicos{Vnt/2)+C2sm{Vnt/2)) = Ae'*^^ cos{\/nt/2- (p). 

Since the roots have nonzero imaginary part, the system is underdamped. 

The damped angular frequency is a;^^ = \/TT/ 2. 

The initial conditions are satisfied when ci = 1 and C2 = 1/ y/li. So, 



x{t) = e-^'^ ( cos( Vll il2) + — ^ sin(Vll tl2) 

Vll 



where ^ = tan-i(l/\/Tl) 



^^e-*/2^os(VlH/2 



2 
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Figure 1: The damped oscillation for example 1. 
Case (ii) Overdamping (distinct real roots) 

If > Amk then the term under the square root is positive and the char- 
acteristic roots are real and distinct. In order for > 4mk the damping 
constant b must be relatively large. 

One extremely important thing to notice is that in this case the roots 
are both negative. You can see this by looking at the formula (2). The 
term under the square root is positive by assumption, so the roots are real. 
Since b^ — Amk < b^ the square root is less than b and therefore the root 
—b + \/b'^ — 4mk < 0. The other root is clearly negative. 

Now we use the roots to solve equation (1) in this case. 

^, -b + Vb'^ - 4mk -b - \/b^ - ^mk 

Characteristic roots: ri = , r2 = . 

2m 2m 

Exponential solutions: e''^^ e''^^ 
General solution: 

Let's analyze this physically. When the damping is large the frictional 
force is so great that the system can't oscillate. It might sound odd, but 
an unforced overdamped harmonic oscillator does not oscillate. Since both 
exponents are negative every solution in this case goes asymptotically to 
the equilibrium x = 0. 

At the top of many doors is a spring to make them shut automatically. 
The spring is damped to control the rate at which the door closes. If the 
damper is strong enough, so that the spring is overdamped, then the door 
just settles back to the equilibrium position (i.e. the closed position) with- 
out oscillating -which is usually what is wanted in this case. 

Example 2. Show that the system x + Ax + 3x = is overdamped and 
graph the solution with initial conditions x{0) = 1, i:(0) = 0. Which root 
controls how fast the solution returns to equilibrium? 
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Solution. Characteristic equation: + 4s + 3 = 0. 
Characteristic roots: (this factors) —1, —3. 
Exponential solutions: e~*, e~^K 
General solution: 

x{t) = cie'^ + C2e~^K 

Because the roots are real and different, the system is overdamped. 

The intial conditions are satisfied when ci = 3/2, C2 = —1/2. So, x{t) = 

3e-72 - e-372. 




01234 
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Figure 2: The overdamped graph for example 2. 

Because goes to more slowly than e~^^^^ it controls the rate at which 
X goes to 0. (Remember, it is the term that goes to zero slowest term that 
controls the rate.) 

Case (iii) Critical Damping (repeated real roots) 

If = Amk then the term under the square root is and the characteristic 
pol5m.omial has repeated roots, —b/2m, —b/2m. 

Now we use the roots to solve equation (1) in this case. We have only 
one exponential solution, so we need to multiply it by t to get the second 
solution. 

Basic solutions: e'^*/^"", te'^*/^"". 
General solution: 

x{t)=e-^^/^"'{ci + C2t). 

As in the overdamped case, this does not oscillate. It is worth noting 
that for a fixed m and k, choosing b to be the critical damping value gives 
the fastest return of the system to its equilibrium position. In engineering 
design this is often a desirable property. This can be seen by considering 
the roots, but we will not go through the algebra that shows this. (See figure 
(4).) 



4 



Under, Over and Critical Damping 



OCW 18.03SC 



Example 3. Show that the system x + Ax + Ax = 
graph the solution with initial conditions x{0) = 

Solution. Characteristic equation: + 4s + 4 
Characteristic roots: (this factors) —2, —2. 
Exponential solutions: (only one) e •^^ 
General solution: 

x{t) = e~^^{ci + C2t). 

Because the roots are repeated, the system is critically damped. 
The intial conditions are satisfied when Ci = 1, C2 = 2. So, 

x{t) =e-2*(l+20. 
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Figure 3: The critically damped graph for example 3. 

Notice that qualitatively the graphs for the overdamped and critically 
damped cases are similar. 

The following figure shows plots for solutions to x + bx + x = with 
initial conditions x(0) = 1, x(0) = 0. The three plots are b = 1 under- 
damped; b = 2 critically damped (dashed line); b = 3 overdamped. Notice 
that the critically damped curve has the fastest decay. 
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Figure 4: Plots of solutions tox + bx + x = 0. 



is critically damped and 
= 1, i(0) = 0. 

= 0. 
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